Fast entanglement of two charge-phase qubits 
through nonadiabatic coupling to a large junction 
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We propose a theoretical protocol for quantum logic gates between two Josephson junction charge- 
phase qubits through the control of their coupling to a large junction. In the low excitation limit 
of the large junction when Ej S> E c , it behaves effectively as a quantum data-bus mode of a har- 
monic oscillator. Our protocol is efficient and fast. In addition, it does not require the data-bus to 
stay adiabatically in its ground state, as such it can be implemented over a wide parameter regime 
independent of the data-bus quantum state. 



PACS numbers: 85.25.Dq,03.67.Lx, 03.65.-w 

Significant progress has been made in improving the 
quantum coherence of Josephson junction (JJ) based 
qubits [1,2] since the first experiment breakthrough in 
1999 [3]. The coherent interaction of two J J qubits has 
been implemented and the initial indications of their en- 
tanglement have been detected [4,5]. Even the demon- 
stration of 2-bit conditional gate of charge qubits was re- 
ported most recently [6] . These developments have paved 
the way towards the realization of the two important el- 
ements of universal quantum computation; the ability 
to implement arbitrary single-bit rotations and the con- 
trolled logic gates between two qubits. Various theoret- 
ical schemes have been proposed for quantum logic gate 
operations of JJ qubits. The most important measure of 
their success depends on the proposed mechanism for a 
controlled coupling between the JJ qubits. Most proto- 
cols involve the coupling of each JJ qubit to an auxiliary 
data-bus (such as variable transformer) and affect an ef- 
fective coupling between two qubits through the elimina- 
tion of dynamic variables of the data-bus [7-13]. 

The existing schemes [7-13] usually apply to weak cou- 
pling parameters, thus leading to slow gates with long 
operation times. Most require the data-bus to be adia- 
batic, i.e. to stay in its ground state or a different pure 
state during gate operation. In this article, we suggest 
a protocol capable of efficient and fast gate operations 
between two charge-phase qubits [2]. As will be shown 
in detail below, our protocol is insensitive to the state of 
the data-bus, thus is not restricted to the weak coupling 
limit as required by the adiabatic condition. In fact, the 
resulting strong effective coupling leads to short gate op- 
eration times. Furthermore, all control parameters in our 
protocol can be modulated within current experiments. 

Our main idea is to construct a system of two J J qubits 
linked by a quantum data-bus of a harmonic oscillator 
with the creation (annihilation) operator a (an. Each 
of the two qubits is linearly coupled to this data-bus 
mode and described by the following form of the Hamil- 
tonian H (t) = aF(t) + h.c. with a general force F(t) = 
(3\{t)Fi + P2(t)F 2 . Fj (j = 1,2) is a dynamic variable 



(not necessarily an observable) of the j-th qubit. It is 
known that when Fj commutes with its own conjugate 
Fj-, the time evolution of H(t) is given by 



U(t) = exp 



i^Vkji^F, 



e a{t)a e° 



(1) 



where the expressions for Hkj{t) and a(t) can be explic- 
itly obtained using the Wei-Norman algebra [14,15]. The 
crucial observation is, at certain instants of time t = T, 
a(T) = and thus the time evolution is simply described 
by 



U(T) = exp 



(2) 



Formally, the quadratic terms of Fj correspond to the ef- 
fective nonlinear interactions between the two qubits. It 
is important to note that even for a fixed time T, there 
still exists enough parameters in the system to modu- 
late Hkj{T) to affect the required logic gate. Clearly our 
strategy does not require the state of the data-bus to be 
adiabatic. Similar approaches have been adopted before 
in protocols of "quantum computing with thermal ions" 
[16] and by Wang et. al. [17]. 
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FIG. 1. Schematic diagram of Two charge-phase qubits are 
coupled through a large junction. 
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We consider a system of coupled charge-phase qubits 
as illustrated by the electronic circuit shown in Fig. 1. 
Similar systems have been recently discussed by You et. 
al. [12]. Our idea works in the limit when the large junc- 
tion (of capacitance C) stays at a low excitation or even 
a thermal state. The same model with only one qubit 
was used to demonstrate the progressive decoherence by 
us [18]. Each of the Cooper pair box is split into two 
small junctions of capacitance C k and C k {k = 1,2) that 
form a superconducting loop. C gk is the capacitance of 
the gate, and Ej k is the Josephson coupling energy. For 

simplicity we set C k — C k = C k and E Jk = E Jk = Ej k . 
If C is much larger than all other capacitances in the cir- 
cuit, then the Coulomb energy can be approximated by 
a non-entangled form. The total Hamiltonian containing 
the reduced Coulomb energy and the three Josephson 
coupling energies, then reads 



where 



H= ^2 [ E ck(n k -n gk f - E Jk (cos ip' k + cos if k ^j 



fc=l,2 



+E C N 2 - Ej cos ( 



(3) 



where E ck = 2e 2 /C sk , n gk = C gk V gk /2e, E c = 2e 2 /C 1 
and Csfc = C gk + C k + C k . Here n k is the number of 
Cooper pairs on the k-th island, while N is the number 
of Cooper pairs on the Coulomb island connected with 
the large junction. ip kl ip kl and 9 superconduction phase 
differences across the relevant junctions. They are related 
through the fluxoid quantization condition around the 
loop 9 + ip k —tp k = 2Q k and Q k = n^ xk /4>o- Introducing 
ipk = (<Pfc + ¥?/J/2, we rewrite Eq. (3) to 



H = J2 E ck(n k - n gk f + E C N 2 - Ej 

fe=l,2 

-2 E Jk cos(--O k ) cos ip k , 

fc=l,2 



cosf 



(4) 



now with the quantization condition [ip k , m] — iSi k . As is 
well known for charge-phase qubit, when n gk — 0.5 and 
when Ej k is not much larger than E ck , the linear com- 
binations of the two lowest charge eigenstates {|0)j,|l)j} 
for each of the split Cooper pair box consist a good rep- 
resentation of a qubit. 

We now consider the "coherent" regime when Ej 3> 
E c . As was found before in the study of quantum phase 
transitions [20] and also some other references such as 
[12], the spectrum of the low energy part of the large 
junction can be described approximately by a harmonic 
oscillator. Within this approximation, we expand the 
above Hamiltonian around 6 = up to 0(9 2 ) and obtain 
an effective spin-boson Hamiltonian 

H = H(g k ,f k ) 

= [9k(a) + a)a xk + f k a xk ] +VLa)a, (5) 

k=l,2 



fl = \J 2E c Ej 
E, 

9k 



32Ej 
fk = -E Jk cos0 fc , 



E Jk sinOfe, 



(6) 



and the quasi-spin and boson operators are defined as 



cr xk 



l>*fc<0| + |0) fcfc <l 
8E r 



e + i 



2E> 



N. 



(7) 



The validity of this approximation can be verified 
through a straightforward numerical calculation. In fact, 
we first calculate the evolution wave function \tpi) , which 
is governed by the approximate Hamiltonian (5). Then 
we compare it with the exact wave function \i[>2 } gov- 
erned by the exact Hamiltonian (3) to a sufficiently high 
precision in expansion of the cosine function of 9 (eg, up 
to 0(9 6 )). In this way we showed that, during one gate 
operation time ( ~ 10~ 10 s as calculated in this paper 
later), the deviation of |V>i) from |-0 2 ) is 10~ 4 at most. 

This equivalent spin-boson system can be solved ex- 
actly since the interaction term g k (a) + a)<r xk commutes 
with the free spin-part. For each spin eigenstate, the a xk 
acts as a linear external force on the boson part. With 
the Wei- Norman algebraic method [14], we find in the 
interaction picture 

Ui{t) = Ui[t\g k \ = exp[-iC(i) - iA(t)a xl <r x2 ] 

x II ^P[-i B k(i)aa x k]exp[-iBk(t)a^a x k], (8) 

fc=l,2 

where the time-dependent parameters are given by 



B k {t) 



9k t-iSU 



-in 



(e- mt - 1) , 



C(t) 



n 

si +52/1 



(9) 



When B k (t) = 0, the evolution operator Ui (t) = 
exp[— iC (t) — iA(t)a x \<j x i\ becomes independent of vari- 
ables of the large junction. It takes the canonical form ca- 
pable of two-qubit gate operations. We note that Bi(t) is 
a periodic function of time and it vanishes at t n = 2mr/tt 
for integer n = 0, ±1, ±2, • ■ -. At these instants of time, 
the time evolution operator becomes explicitly as 



Ui (i„) = cxp 



Anirgig 2 



-O x \O x 2 



up to a phase factor exp [—i2nir(g 2 + g 2 )/i} 2 ] in the in- 
teraction picture. This is equivalent to a system of two 
coupled qubits with an interaction of the form cx a x 1 er^ • 
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For a two-qubit system governed by a Hamiltonian 
of the form g<J x \o x2 , one usually adjusts the evolution 
time to realize an arbitrary two-bit controlled rotation 
U X x (£) = exp[-^cr x iCT x2 ] at t = £,/g. In our system, 
the evolution time is fixed at t n by the requirement of 
Bi{t n ) = 0. This is in fact not a problem as experi- 
mentally one can vary g\ and 172 to affect the desired 
rotation U xx (^) at t n . We note that gi depends on the 
external flux <§> x k of the loop 6^ = n& x k/<f)o, so it is easy 

to adjust to the maximum of gtmax = gk (@fe = f ) = 

v 1 

E, 



MM' 

\32Ej J 



min — 



32Ej 1 Ejk- The minimal time for one operation 
" " n,7r is given by 



n 



n 2 



2 5i(f)32(f) 



+ 1, 



(10) 



where [x] denotes the integer part of x. 

At the same time, we note that the free Hamiltonian 
part Ho = fk&xk does not vanish during the above 
discussed two-qubit operation. However, it commutes 
with the interaction term oc u x \ cr X 2 and simply leads the 
evolution operator being 



U xx (£) = jg-i^i^^g-ffi/fcdV^fcjgiE/fcdJa-a, 



(11) 



i.e. augmented by two single-bit operations with (£) = 
—Ejk cos0fc (£). 9fc (£) satisfies the equation 

1 

mr ( 8E \ 2 

— f-^j Ej 1 Ej 2S mQ lS mQ 2 = Z. (12) 

Alternatively we can remove the influence of the large 
junction on the two-qubit logic operation by using two 
operations in succession, a concept as presented recently 
for a two-qubit gate on trapped ions [19]. The two steps 
are described as follows. 

First, we evolve the system with the Hamiltonian H 
for a duration r/2. The evolution operator (in the inter- 
action picture) is 



(3) 



Ui 



2' 9k 



(13) 



Second, at time t = r/2, wc instantly reverse the direc- 
tion of the magnetic field such that a sudden change of 
flux from <& x k to — <$> x k occurs that leads to g k = —gk and 
f k = fk ■ The system is then driven by a new Hamiltonian 
H' = H (—gk, fk), with which we evolve for another r/2. 
Since Hq = H and Hj = —H T , the evolution operator 
(in the interaction picture) becomes 



U 



Ui 



'9k 



(14) 



The combined dynamics from the above two steps is 
now described by the following evolution operator 



= Ui (h- gk ) Ui (h 9k ) 

= exp[-iM(T)(T xl a x2 ] exp[-iN ( T )] ; 

where 



(15) 



M{t)=aC-)+A( T -) 



2glg2 

n 



2 n T 



(16) 



and 



N(t)=d'C-) + D( T -) 



9l 



g\ ( 2 . nr 
" ' — sin — 

n 2 



(17) 



Again we see that after the above two successive oper- 
ations, we realize an effective controlled rotation of the 
two qubits, i.e. U xx (£) = exp[— i^a x \a x2 \ can be real- 
ized by fixing Qj~ = ir/2 or Q x k = 4>q/2 and for a time 
t satisfying M(r) — S,. As a bonus we find that fk = 
in this case, i.e. the single-bit rotation vanishes auto- 
matically during the two-bit rotation. With both of the 
above approaches, single bit operations can be cleanly 
implemented by setting Qk = ^ or 9k = 0- 

We now consider the implementation of our protocol in 
realistic experiments. We note the typical time for gate 
operation is 



27rf7 



2tt 

— = P — > 
9 9 



(18) 



with g to be roughly understood as the order of mag- 
nitude of gk and 2ir/g approximately the single qubit 
operation time. One can increase the ratio p — g/Vl to 
shorten these operation times similar to other JJ cou- 
pling schemes [7 12]. It is important to emphasize, how- 
ever, most of these other schemes are based on an adi- 
abatic evolution of the data-bus, thus are limited to a 
weak coupling, or, a small p. In contrast, our protocol is 
not confined to the adiabatic dynamics, thus can operate 
with a larger p for a faster gate. In fact, even when the 
capacitance of the large junction is much larger than the 
other capacitances, the ratio ^ can still be made large 
with realistic experimental parameters. For example, if 
wc take E.j = 800 (/zeV), E c = 10 (pcV), E ci = 200 
( M cV), Ejk = 200 OueV), $ xfc = o /2 [2,21], wc find 
g/fl ~ 0.23, a limit may be prohibited for other schemes 
yet works well for our protocol. 

A key requirement for our protocol is the precise con- 
trol of the operation time, especially in the second ap- 
proach which needs a sudden switch of the magnetic field 
during each two-qubit operation. However, the second 
scheme need not to change the magnitude of the mag- 
netic field but just to switch the direction of it. Before 
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concluding, we wish to remark that, although our circuit 
resembles that of You et. al. [12], our operating scheme 
and the underlying physics is different. The essential dif- 
ference is that our gate does not require an adiabatic 
operation, thus is faster. 

In conclusion, we have presented an efficient proto- 
col for implementing controlled interactions between two 
charge-phase qubits. We have discussed two alternative 
approaches to realize our protocol. Our scheme seems 
to work over a wide parameter range and is faster than 
most existing protocols. From the experimental point of 
view, our protocol seems advantageous as it only needs 
the control of one system parameter, that is the external 
applied magnetic field. We have adopted a setup involv- 
ing charge-phase qubit, thus our protocol is less sensitive 
to charge fluctuations along with phase fluctuation [2]. 
In addition, our circuit arrangement can be easily scaled 
up to larger number of JJ qubits with similar controls. 
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